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The field of plasmonics has experience a renaissance in recent years by providing a large variety of new 
physical effects and applications. Surface plasmon polaritons, i.e. the collective electron oscillations at the 
interface of a metal/semiconductor and a dielectric, may bridge the gap between electronic and photonic 
devices, provided a fast switching mechanism is identified. Here, we demonstrate a surface 
plasmon-polariton diode (SPPD) an optoelectronic switch that can operate at exceedingly large signal 
modulation rates. The SPPD uses heavily doped p-ii junction where surface plasmon polaritons propagate at 
the interface between n and p-type GaAs and can be switched by an external voltage. The devices can operate 
at transmission modulation higher than 98% and depending on the doping and applied voltage can achieve 
switching rates of up to 1 THz. The proposed switch is compatible with the current semiconductor 
fabrication techniques and could lead to nanoscale semiconductor-based optoelectronics. 

The demands for faster data transfer and processing has increased dramatically over the past five decades. This 
need led to enormous advances in the semiconductor industry and a continuous progression towards 
smaller, faster, and more efficient electronic devices \ The scaling of the electronic devices also brought 
about a myriad of challenges. The most daunting problems inhibiting significant increases in the processor speed 
are thermal and resistor- capacitor (RC) delay time, associated with the electronic interconnections and scaling of 
the devices^. Over the past few years saturation in the microprocessor clock speed has been observed due to the 
scaling issues and excessive heat dissipation \ 

On the other hand, optical interconnects possess an unimaginable data carrying capacity, which may allow 
photonic components to mitigate the present bottlenecking in increasing the computational speed^. However, 
their implementation is hindered by the physical size or dimensional mismatch between electronic and dielectric 
photonic components. Moreover it is expected that CMOS foundries will further decrease the feature sizes on 
silicon chips, ultimately down to 10 nm. Such reduction in size would further increase the dimensional mismatch 
between the electronic and optical elements, as the size of the dielectric photonic devices is restricted by the 
diffraction limit^. Moreover, miniaturization introduces several problems including dielectric break down, hot 
carriers and short channel effects, that degrade device reliability^. 

The Surface Plasmon Polaritons (SPPs), spatially confined transverse magnetic (TM) electromagnetic modes 
propagating at the metal -dielectric interfaces, offer the bandwidths of photonic devices and physical dimensions 
shared with nanoscale electronics^'^"^. The potential of plasmonics to bridge the gap between electronics and 
photonics is now well recognized by the scientific community with a large number of investigators working in the 
field of plasmonics^. In 2008 Brongersma et at. demonstrated all-optical switch based on SPP waveguide that uses 
metallic (passive) nanostructures coupled with active PMMA films with photochromic molecules^. 
Unfortunately, the switching rates of the photochromic molecules were low (—20 ns)^. The same group also 
proposed extremely compact gain-assisted plasmonic switch consisting of a gold-air-gold plasmonic waveguide 
side-coupled to a cavity filled with a semiconductor (InGaAsP) gain materiaP^. A metal-oxide-Si field effect 
plasmonic modulators and all-optical modulation by plasmonic excitation of CdSe quantum dots have been 
investigated showing moderate transmission modulation at visible and telecommunication frequencies^ 
Recently, a THz all-optical switch based on a carbon nanotube metamaterial has been proposed, however the 
device shows a rather low transmission modulation of less than 10%^^. 

In this article we propose a new optoelectronic device which we refer to as a Surface Plasmon Polariton Diode 
(SPPD). The SPPD demonstrates active control of SPPs at the interface between doped n-type and p-type 
semiconductors. While the propagation of SPPs on metal surfaces is well understood, not much attention has 
been paid on the fact that due to the extremely high plasma frequencies inherent to noble metals (for silver cOp = 
9 eV) there is a substantial bandwidths mismatch with electronics that cannot be easily bridged. Here we show 
that use of highly doped semiconductor material (GaAs) can serve three distinct purposes: (/) act as a metal-like 
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(SPPD). 



Basic schematic of a Surface Plasmon Polarition Diode 



interface allowing SPP propagation, (//) provide faster electronic SPP 
switching facilitated by a p-n junction, and (///) provides tunable THz 
operational frequency range. More importantly, due to strong local- 
ization of the SPP at the p-n junction, the SPPD can exhibit extremely 
fast (up to 1 THz) switching rates, has relatively small size and can be 
used to develop the fundamental logic for prospective use in optoe- 
lectronic data processing. 

Figure 1 shows a basic schematic of the proposed SPPD con- 
structed of a highly doped p-n junction formed between two control 
electrodes. The SPPD basic principle of operation follows. For fre- 
quencies below the plasma frequency, cOp, the semiconductors 
behave as metals and can be tuned by doping or by applying an 
external bias. The metal-like properties of the semiconductors allow 
propagation and active control of SPPs at the p-n junction. Under 



zero-bias the SPPs propagate freely at the metallurgic junction and 
across the drift- diffusion zone, estabUshing the ON-state of the 
devices. With the introduction of the external bias electrons are 
injected in the p-type region, dramatically altering its optical prop- 
erties from dielectric-like to metal-like. This prohibits SPP propaga- 
tion across the active zone and establishes the OFF-state of the 
device. The SPPD can be designed using a large variety of semi- 
conductor materials. Here we use GaAs instead of Si due to its low 
effective electron mass, higher plasma frequency, and high mobility. 

Figure 2(a) depicts the SPPs dispersion curves w = w(kspp) (where 
kspp is the SPP wavevector) for two doping concentrations (see 
Methods). At low frequencies the SPP wave vector coincides with that 
in air while at frequencies close to the surface plasmon frequency 
cOsp = cOp / \^28h {sb — 12.9 is the bound electrons contribution to 
the p-type layer permittivity^^) it experiences a "resonant" type of 
behavior with kspp:>2K/X, where X is the wavelength in free space. 
Consequently, the effective mode wavelength is dramatically reduced 
below the diffraction limit in the semiconductor, allowing for the 
SPPD to be substantially smaller compared to similar devices based 
on bulk optical modes. We must note that when choosing the SPPD 
operation frequency it is important to consider the tradeoff between 
device size and mode dissipation. In this work we frx the operation 
frequency at 35 THz which allows for a substantial decrease of the 
mode size less than 300 nm in a direction perpendicular to the p-n 
junction, without imparting extensive propagation loses. Further- 
more, at low doping concentrations the depletion region can act as 
a gap waveguide for weakly bounded bulk modes^^. For the doping 
levels considered in this work (Nd > 10^^ cm~^) the depletion region 
is only a few nanometers in size^^ and the excitation of such weakly 
bound photonic modes is strongly inhibited (see Fig. 2(a)). 

The SPPD operation characteristics are first studied under steady 
state conditions. We use a self-consistent numerical model that 
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Figure 2 | SPPD steady-state response, (a), The SPPD dispersion curves in air/p-type/n-type GaAs waveguide configuration for two doping 
concentrations; Na = 1 X 10'\ = 5 X 10'' cm~^ (red line), and Na = 1 X lO''', = 3 X 10'^ cm~' (blue line), (b), Steady state minority carrier 
concentration profiles for different applied voltages, (c), The transmittance curves obtained using the integrated finite difference semiconductor and 
optical module (solid blue line) and the WKB approximation (dashed red line), show exponential decrease of the output signal for V> = 1.57 V and 
(d), SPP propagation across the SPPD drift- diffusion zone obtained at different applied voltages. In the calculations the operation frequency is set at 
35 THz, the thickness of the p-type layer is d = 0.5 fim and the overall length of the active drift- diffusion region is fixed at w = 2 jum. 
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Figure 3 | SPPD switching, (a), SPPD transmittance curve (solid blue line) under step-type of input voltage (dashed red line). The signal is repetitively 
switched following the external voltage, (b), Electron concentration and (c), Magnetic field across the device at two different times corresponding to the 
OFF and ON states of the device. In the calculations the operation frequency is set at 35 THz, and the doping concentrations are = = 5 X 

10^'^ cm-\ 



merges the electronic response, obtained from finite difference (FD) 
integrated circuits semiconductor code (SENTAURUS, Synopsys 
Inc.), and the optical response of the device simulated with a com- 
mercial electromagnetic code (COMSOL Multiphysics). Under zero 
apphed bias the electron concentration in the p-type layer is low 
n—no = n] /Na, where Ui = 2.25 X 10^ cm~^ and Na is the acceptor 
concentration, causing this layer to behave as a dielectric (See 
Fig. 2(b) )^^. Concurrently, the uninterrupted SPP propagation estab- 
lish the "ON" state of the device (See Fig. 2(c,d)). If forward bias is 
applied the electrons drift from the n-type layer into thep-type layer, 
and exponentially increase the minority concentration in close prox- 
imity to the metallurgic junction /i~/io^^^^^, where Vt = ksT/q = 
0.026 V is the thermal voltage). In the presence of sufficient bias 
{V > Vc) the electron concentration in the drift- diffusion zone can 
surpass a critical value at which point the permittivity of the/>-doped 
layer becomes negative (Spico) < 0), and this layer becomes a metal- 
like (the semiconductor permittivity is described by the Drude 
model, see Methods). The SPP could no-longer propagate through 
the active drift- diffusion zone, and the OFF state of the device is 
established (Fig. 2(c,d)). For a given operation frequency, w, the 
critical voltage can be obtained from the transparency condition 
w = Wp{Vc) / \f&h (or ep(<^j V^c) — 0)» ^nd the carrier balance equation 
n = fic^fiQe^'f^^ , giving 



where co^o ^ (^o^^/^o^e)^^^ is the plasma frequency of the minority 
carriers under thermal equilibrium and zero external bias, is the 
electron effective mass in the p-type layer, and q is the electron 
charge. For AT^ = 5 X 10^^ cm~^ and/= ojlln = 35 THz the critical 
voltage is ^ 1-57 V, which is more than an order of magnitude 
lower than the breakdown voltage for the device^^. Due to the low 
effective mass of electrons in the/>-doped GaAs (m^ = 0.067 mo) the 
critical voltage is lower than expected for Si-based (m^ = 0.26 mo) 
SPPDs, which leads to lower heat dissipation and further justifies the 
choice of using GaAs. 

The SPPD switching, can be qualitatively described using the 
Wentzel-Kramers-Brillouin (WKB) method^^, which implemented 
for the SPP transmittance across the drift- diffusion region yields: 

T{V) = Jo"^^^^'('''^)'^''|^^e"^'^^"'^^^^^^l (2) 

where w is the length of the active drift-diffusion region (roughly 
corresponding to the length of top electrode), and the SPP wave 
vector, kspp, which is generally position and voltage dependent due 
to the spatially inhomogeneous minority carrier concentration in the 
p-type layer. Under steady state conditions, the SPP wave vector in 
the drift-diffusion zone can be assumed to be spatially homogeneous 



and dependent only on the applied voltage kspp{Xy V) ~ KppiV). As a 
result, the WKB approximation gives an explicit, analytical formula- 
tion of the SPP transmission across the device. For a given frequency 
of operation, w, and sufficiently low minority carrier concentration 
ksppEU and T ^ 1. When V > V^, the permittivity of the p-doped 
layer becomes negative (Spiw) < 0) and the SPP wave vector is 

complex quantity (kspp'^^'i ^^^^^ e^eC). Consequently the transmit- 
tance exponentially decreases with an increase in the applied voltage 
(see dashed line in Fig. 2c). 

The prime interest in the study of the SPPD is the speed at which 
the optical properties of the/>-layer switches between dielectric (Sp > 
0) to metallic {Sp < 0) and vice-versa. The switching rate depends on 
the applied bias and doping concentrations and can be obtained from 
time-dependent calculations of the minority concentration and the 
corresponding SPPD propagation characteristics. Figure 3(a) shows 
the SPPD temporal responses under step-type input voltage where a 
forward bias is applied across the drift- diffusion zone for a period of 
5 ps and a zero bias follows for 10 ps (see Fig. 3(a)). The SPPD 
transmittance curve, depicted in Figure 3(a), is obtained by FD 
self- consistent numerical calculations in the time domain and shows 
reproducible switching. Under a forward bias larger than the critical 
{V> Vc), the minority carries easily overcome the potential barrier 
due to the space-charge region. The electron concentration rises 
exponentially close to the metallurgic junction. Within a few pico- 
seconds the concentration surpasses the critical value, n^, at which 
point the SPPs are no longer able to propagate (establishing the 
"OFF" state of the device), see Fig. 3(b)&:(c). Once the applied bias 
is removed the excess electrons in the />-layer diffuse and their con- 
centration falls below the critical value in approximately 5 ps. At this 
point the SPPs can propagate through the device re-establishing the 
"ON" state. The results indicate that extremely fast switching can be 
achieved corresponding to data rates of 100 GHz and faster. At the 
operational frequency of 35 THz and a voltage of y = 4 V the SPPD 
transmission modulation is moddB =10 logio(^'majc/^m/n) = 18 dB 
(corresponding to 98.4% transmission modulation), surpassing by 
far what have been reported in the literature, in the case of all-optical 
SPP based switches^'^'-^l 

From Figure 2(b) and Figure 3 (b) it is evident that it is not 
necessary for the minority concentration to surpass the critical value 
within the entire drift- diffusion zone to switch the SPPD, rather, only 
in close proximity to the metallurgic junction. This is a fundament- 
ally different effect compared to conventional, all -electronic diodes, 
where the response time is proportional to the electron fly time 
between the terminals. The SPPs, however, are exponentially sens- 
itive to the dielectric environment in close proximity to the space- 
charge region, the metallurgic interface, and are switched OFF within 
a few picoseconds once the electron concentration within 20 to 
50 nm from the junction surpasses the critical concentration value. 
He. This fact and the almost immediate signal transfer between a set of 
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SPPDs (the SPPs travel with the speed of hght) is what may allow for 
much faster response times and negligible interconnect delay times 
of lumped optoelectronic circuits made of SPPD as compared to 
conventional diodes. 

The temporal response of the device depends on the injection rate 
of electrons from the n-type into the p-type layer and thus on both 
applied voltage and doping concentrations. Figure 4(a) and 4(b) 
present parametric studies of the SPPD response times as function 
of forward bias and p-doping concentrations. The ON time of the 
device is diffusion limited and is thus independent of the initial 
applied bias. The OFF times, however, are governed by the drift of 
the electrons in the p-type layer and are thus inversely proportional 
to the applied voltage. We must caution that operation at high 
external biases is not encouraged since it will lead to substantial 
increase in power dissipation and potential device failure. Thermal 
analysis of the SPPD has to be performed to establish the maximum 
voltage that can be applied, which will be the subject of future work. 
Furthermore, as shown in Figure 4(b), the response times decrease 
exponentially with increase in thep-layer doping. This phenomenon 
is due to the dependence of Tqn and Tqff on the SPP penetration 
depth in thep-type layer which varies with the doping. The effects of 
applied voltage and doping on the response times can be quantita- 
tively estimated as follows. 

Due to the presence of external bias, electrons drift into thep-layer 
with a velocity //^{V — Vi,i) / Xp, where f/^ is the minority carrier 
drift mobility in/>-type GaAs^^, Xp is the thickness of thep-layer and 
Vi,i is the build-in potentiaPl The SPPD OFF time is inversely 
proportional to the drift velocity and the applied bias and can be 
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estimated as Tqff = Uppl^d^ where kpp = l//cf ^ = 1^ y av-^^ - tpAVg 

the field penetration depth of the SPP in the p-type layer (see 
Methods). At zero bias, the excess minority carriers in the p-layer 
are removed by the net carrier outflow from the quasi-neutral region 
via diffusion. Using dimensional analysis we can write the response 

(ON) time as tqn = l^ppj {^FDn), where F ~ 4 is a shape factor and it 

is near-independent on the parameters under investigation, 
Dn = {kT/q)fi^ is the diffusion co-efficient of the electron in p- 
region, fi^ is the electron mobility^". The predicted response times 
by the above analytical considerations are shown with solid lines in 
Figure 4 and can be considered as the lower bounds of the actual 
response times. Overall, the SPPD can operate at very fast switching 
rates of up to 1 THz (or 1 ps response time). 

The SPPD, similarly to its all -electronic counterparts, can be 
potentially used to develop the fundamental unit of a microprocessor 
i.e., a logic gate. Due to the fast response times the SPPD may sig- 
nificantly improve the clock speed of an optoelectronic processor 
unit. In Figure 5 we propose and study possible designs of optoelec- 
tronic logic elements based on the SPPD. The logic elements are 
constructed by connecting three SPPs waveguides in a "Y" shape 
configuration in which each waveguide is controlled by a separate 
SPPD. The NAND gate (Figure 5(a)) is controlled by two external 
voltages, Vi and V2, applied at the input waveguides SPPDs, while the 
NOR gate (Figure 5(b)) is controlled by external voltages applied in 
series at the input and output waveguides SPPDs. In the absence of 
external bias {Vi = 0 V and V2 = 0 V) the input SPPs interact 
constructively at the Y- junction, and the signal is transferred estab- 
Ushing the "ON" states for both NAND and NOR gates. In the 
presence of single external bias (Vi = 0 y, ^2 — 2 y or Vi = 2 y, 
V2 — 0 the SPPs from one of the impute waveguides is transmit- 
ted in the case of the NAND gate ("ON" state) and reflected in the 
case of a NOR gate ("OFF" state). When both input SPPDs are 
biased, there is no transmission and both gates are in "OFF" state. 
The actual three dimensional magnetic field profiles calculated in 
correspondence to the logic tables Fig. 5 (c, f) are depicted in 
Fig. 5(b, e). Both electronic control (through the SPPD applied vol- 
tages) and constructive/destructive interference at the output wave- 
guide are utilized to establish the proper input/output signal 
associations. One limitation for the proposed "Y" shaped device is 
it's the phase sensitivity. Specifically, the gate requires two input 
signal of same phase to have a constructive interference in the 
absence of external bias (Vi = 0 V and V2 = 0 V) and thus the 
highest power output. We must emphasize that in principle there 
are many possible gate configurations and the chosen "Y" shape 
design may not be the optimal. More importantly, due to the 
SPPD nature of operation, the proposed logic gates are optoelectro- 
nic (not all-optic) and can be simply viewed as a fast and efficient way 
of controlling the propagation of surface plasmon polariton (SPP) 
modes at highly doped semiconductor interfaces. 

In conclusion we have introduced and demonstrated a new optoe- 
lectronic switching device, i.e. surface plasmon polariton diode 
(SPPD), capable of operating at exceedingly high switching rates 
and signal modulation. The fast device response is achieved by active 
control of the propagation of SPPs along a highly doped p-n junction. 
The performed steady-state and time dependent numerical analysis 
suggests that the SPPD can operate at signal modulation in excess of 
18 dB and switching rates from 100 GHz up to 1 THz. Particular 
designs of the fundamental logic elements ("Y" shape configuration), 
NAND and NOR gates, have been demonstrated. 

Methods 

Surface plasmon polaritions in doped semiconductors. Similar to metals, the 
optical response of highly doped GaAs can be characterized by the Drude model with 

permittivity given as sgms (c^) = ^b — (^) / + i(i>cOj{n)), where = 12.9 is the 

bound electron permittivity^^. An important difference to note is that both the plasma 
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Figure 4 | SPPD response times, (a), Dependence of the SPPD ON and 
OFF times on the applied voltage {V> V^). The device OFF time is 
dependent on the electron drift toward the electrodes and is thus inversely 
proportional to the applied voltage. In the calculations the SPP operation 
frequency is set at 35 THz, and the doping concentrations are fixed at = 
No = 5 X 10^*^ cm~^. (b), The SPPD response times as function of the p- 
layer doping concentration. The doping influences both the width of the 
space charge region and the electron mobility thus modifies to a different 
extent the ON and OFF times of the device. In the figures the numerical 
data (diamonds) is compared to the analytical drift-diffusion model (solid 
lines). In the calculations the applied voltage is ^ = 3 V. 
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frequency, cOp = (ne^/EQmeY'^, and relaxation rate, co-, = e/mefi{n), now depend on the 
doping concentration. This dependence allows for fine tuning of the optical 
properties of the doped GaAs for frequencies up to 200 THz (corresponding to a 
maximum doping equal to 6 X 10^^ cm~^ ^^). The SPPs dispersion relation^^ for a 
three layer planar waveguide configuration (see Fig. 1) is given as 

^-4k,d_ kd/£d-\-km/Sm kd/Sd-l-kg/Sg 
kd/^d - kml&m kdl&d-kal&a 

where e^^, are the dielectric permittivities of ^-layer, p-layer and surrounding 

dielectric (air), d is the thickness of thep-layer, ki = ,J k'^^p — sikl, lE{m,d,a}. The 

SPPs wave vector k^pp is larger than that in free space ko = (co/c) (see Fig. 2(a)). 
Therefore, coupling and de-coupling of energy with the far field is inhibited due to the 
momentum mismatch. Consequently, the SPPs are exponentially confined to the p-n 
junction interface^'^. 

Numerical calculations. The SPPD's operation is been modeled using finite- element, 
steady state and time domain numerical calculations. A self- consistent model was 
developed by first calculating the electron and hole concentrations over the SPPD 
using commercial semiconductor device simulator (SENTAURUS, Synopsys Inc.). 
The inhomogeneous carrier's concentrations for the n-layer and p-layer (including 
the depletion region concentration) are then used to extract the position dependent 
effective permittivity of the n and p-layer. Finite element calculations of the 
electromagnetic response of the diode based on COMSOL Multiphysics, follow 
accordingly. This process is repeated for various input parameters such as donor/ 
acceptor doping concentrations, applied bias, and operation frequency. A common 



(MATLAB) environment is used to merge the two commercial codes, which allows 
automatically sharing data and performing large sets of parametric studies, critical for 
the optimization of the SPPD performance. The response times are estimated by 
fitting the transmittance data (similar to Fig. 3a) with an exponent e~^'\ where t is the 
switching time. 
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